Quantum phase transition in the Hartree-Fock wave function of the hydrogen molecule 
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Precise solutions of the Hartree-Fock equations for the ground state of the hydrogen molecule are 
obtained for a wide range of internuclear distances R by means of a two-dimensional fully numerical 
mesh computational method. The spatial parts of the single-electron wave functions are found to 
be coinciding for R < 2.30 a.u. At larger distances they become different and as R — > oo each of 
them takes the form corresponding to a separate atom. This quantum phase transition occurring at 
R = 2.30 a.u. gives a natural boundary between a delocalized behavior of electrons in the molecule 
and their distribution over separate atoms. This phenomenon can be also considered as an analog 
of the Wigner crystallization or the Mott transition on the level of a single molecule. 



In recent years multiple studies of quantum phase tran- 
sitions were carried out. These studies are significant 
for understanding of many complicated phenomena in 
solid states, clusters and so on. As an example, stud- 
ies of the Wigner crystallization in two-dimensional sys- 
tems can be presented Q. Many of these effects are as- 
sociated with symmetry breakings in the Hartree-Fock 
ground states. On the other hand, spontaneous symme- 
try breakings in molecular systems were reported mainly 
in a negative context H . In this communication we study 
the simplest molecule with more than one electron, i.e. 
the H2 molecule, for arbitrary internuclear distances R. 
This molecule was investigated in a number of works by 
means of many precise methods (see Hj). On the other 
hand, there are no studies of precise solutions of the 
Hartree-Fock equations for the H2 molecule for arbitrary 
internuclear distances. We present such calculations and 
show that at some R = R CI the wave function of the 
ground state undergo a quantum phase transition associ- 
ated with a spontaneous symmetry breaking. This tran- 
sition separates a phase with two fully delocalized elec- 
trons with equal spatial parts of their wave functions and 
a phase with two electrons localized on separate atoms 
(i? > R C i ). This phenomenon is very similar to the form- 
ing of Wigner crystals or Mott transitions in the solid 
state context. Due to this phase transition the solution 
of the Hartree-Fock equations has a reasonable physical 
meaning for arbitrary R values, that does not take place 
in a traditional approach with single-electron wave func- 
tions having the symmetry of the molecule as the whole 
at all the R values Q]. This result allows us to expect, 
that, as well as in the solid state physics, some symme- 
try breakings in more complicated molecules and similar 
systems also could be not computational artifacts j2| but 
have an appropriate physical meaning. 

We solve a system of the Hartree-Fock (HF) equations 
by a fully numerical two-dimensional finite-difference 
(mesh) method || ^, fjj. The equations are presented 
in cylindrical coordinates (p, <f>, z) with the axis z coin- 
ciding with the molecular axis. The point z = cor- 
responds to the center of gravity of the molecule. For 



each single-electron wave function a definite value of the 
magnetic quantum number is supposed (m = for the 
ground state, considered below), so that a numerical so- 
lution is carried out on the plane (p, z). Thus, the wave 
function of the spin singlet ground state has the form 

where ip\ and ^2 are solutions of the Hartree-Fock equa- 
tions for the first and second electrons. The correspond- 
ing energy of the electron system we denote as E c . When 

0i 7^ ^2 this wave function is not an eigenfunction of the 
-2 

spin operator S . In this case it is possible to consider 
also a spin-symmetrised two-determinant wave function 

* Symm = bPl (Pl ,Zl)ih (P2 ,Z 2 )+^2 (Pi ,Zl) Ipl (P2 ,Z 2 )] 

and the corresponding energy El ymm . The latter wave 
function is an eigenfunction of the operator S 2 as ipi ^ 

i>2- 

A traditional approach to solution of the molecular 
Hartree-Fock problem employs fully delocalized molecu- 
lar orbitals. In our case this means if>\ = tp2- This allows 
obtaining reasonable results only for relatively small R 
values. For large distances atomic orbitals strictly local- 
ized near corresponding nuclei are sometimes considered. 
Of course, every combination of such approaches gives 
rise to a problem of intermediate values of the parameter 
R. In our calculations we are free from limitations asso- 
ciated with a choice of basis functions and do not require 
this or that way of behavior from ipi and ip2- They are 
direct solutions of the initial equations. 

The results obtained in our numerical solution of the 
Hartree-Fock equations are presented in the Table | and 
Figures p]|| All the data are given in atomic units. As 
one can see, the solutions have very different characters 
for R < 2.3 and for R > 2.3. Several plots of the spa- 
tial parts of the single-electron wave functions and 
■02 for different R are given in Figure § For R > 2.3 
we have ipi ^ ip2 (functions a and b in the plots), and 
with growing R these solutions transform into wave func- 
tions of separate hydrogen atoms. On the other hand, at 
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R < 2.3 the functions tpi and ip2 fully coincide. This 
relation between ipi and ip2 as well as the existence of 
the critical point R = R CT for solutions of the Hartree- 
Fock equations can be illustrated by the behavior of 
their overlap integral {4 , i\ip2} given in Table || and Fig- 
ure H For R < R CI an exact equality (V'ilV'2} = 1 takes 
place, whereas for R above the critical point the function 
{ipi\ip2) {R) shows a near to an exponential decrease with 
increasing R values. Our calculations in a vicinity of the 
critical point allow to estimate its position as R CI = 2.30 
or, more precisely, R CI = 2.29999. For R < R cr we have 
ipi(p,z) = ipi(p, —z) (i — 1,2), whereas for R > R cx this 
symmetry condition does not take place and only con- 
dition ipi(p, z) = ifj2{p,—z) is preserved. It should be 
emphasized, that the transition from the single-electron 
wave functions, symmetric with respect to the center of 
the molecule, to the non-symmetric ones is an intrinsic 
property of the Hartree-Fock equations for the hydrogen 
molecule. Thus, for the hydrogen molecule at R > R CI 
a spontaneously broken symmetry in the Hartree-Fock 
wave function takes place. A similar effect has been found 
previously for the beryllium atom . From the point of 
view of the numerical techniques the Hartree-Fock iter- 
ations can be started from initial wave functions with 
arbitrary small deviations from exact z — parity. The 
final result of iterations does not depend on the degree 
of this initial asymmetry (this could be simply round-off 
errors) and has the form presented in Figure 

The Hartree-Fock electron energies of the hydrogen 
molecule as well as numerical values of the (V'ilV'2) are 
presented in Table Q. For R < i?, cr the Hartree-Fock elec- 
tron energy is an energy calculated on a one-determinant 
wave function of the whole system E c . Available pre- 
cise HF result for the near to the equilibrium distance 
R = 1.4 (E e = -1.8479152858, @) coincide with our 
one. Our energy for R = can be considered as more 
precise than the well known result for the helium atom 
E c = -2.8616799 §. For R > R cr it is possible to calcu- 
late the spin-symmetrised energy Ef™ 1 ™ 1 also presented 
in the table. This energy lies lower than E c except the 
cases (^iIV^) = and (^ll^) = 1, when both energies 
evidently coincide. 

As a discussion it is expedient to note that the result 
obtained above contradicts to a frequently encountered 
opinion, that even the simplest of molecules, the hydro- 
gen molecule, dissociates incorrectly within the Hartree- 
Fock approximation When we do not impose some 
artificial conditions (i.e. the symmetry of the molecule as 
a whole) on single-electron wave functions, the consistent 
Hartree-Fock approach allows obtaining quite reasonable 
solutions for arbitrary internuclear separations. An ex- 
ception are Van der Waals forces, which cannot be ob- 
tained in the Hartree-Fock approximation. In the limit 
R — > 00 our solution describes two separate atoms (H 
+ H) with electronic wave functions tp\ and ^2 respec- 
tively. On the other hand, when establishing symmetry 




FIG. 1: Single electron Hartree-Fock wave functions of the 
hydrogen molecule for different internuclear distances R. For 
R > Rd ~ 2.3 a.u. a (light-gray lines) and b (black lines) 
are the functions with electron densities concentrated at the 
left and right nuclei respectively. For R < R CI these wave 
functions coincide. 
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TABLE I: Physical parameters of the Hartree-Fock hydrogen 



molecule dependent on 


the inter nuclear 


distances. Atomic 


units. 








R 


E c 


nsymm 




0.0 


-2.861679998 




i 
1 


0.2 


-2.7608175 




i 
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0.4 


—2.5798578 




i 
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-2.396660751 




1 
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0.8 


-2.23074400 




1 
1 


1.0 


—2.085138396 




i 
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1.2 


-1.958359581 




1 
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1.3 


-1.901255170 




1 
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1.4 


-1.847915286 




i 
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1.5 


-1.798039362 




1 
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1.6 


-1.751347985 




1 
1 


1.8 


-1.666513946 




i 
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2.0 


-1.591619847 




1 
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2.2 


-1.525186768 




1 
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-1.494754132 


— 1.4y4 ( 04Z0O 


n 0000007 
U.yyyyyy ( 


2.4 
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2.8 


-1.381746342 
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-1.35055131 
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fl ^71 A QOOQ 

U.o / l4yzzo 


3.5 


-1.29277168 
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U.ooDo ( o4o 
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-1.252894522 
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n O^OQOQflf^O 
U.zOzozoUOy 


5.0 


-1.200485964 


1 OflOT/IQQ^f^ 
— l.ZUz ( 4ooDD 


fl 1 01 A QQQ7Q 

U.lz44ooy (o 


6.0 


-1.166747578 


— 1.10 1 zU4zyo 


U.UOOOoO lO( 


7.0 


-1.14287045 


— i.i4zyooyo 


U.UzozU r Ulo 


8.0 


-1.12500216 
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1 111111 A({ 

— l.llllll^iO 


— 1.1111 14zU 


n nn/i o^i ^.on 
U.UU4yolozU 


10. 


-1.10000006 


-1.10000053 


0.00215182 


11. 


-1.09090911 


-1.09090918 


0.0009244618 


12. 


-1.08333334 


-1.08333336 


0.0003931164 


13. 


-1.07692309 


-1.07692309 


0.0001656485 


14. 


-1.0714286 


-1.0714286 


0.0000692334 


15. 


-1.0666667 


-1.0666667 


0.0000287261 



conditions tpiifl, z) = ijji(p, —z) or/and = ip2i we have 
in the limit R — > oo a traditional non-physical Hartree- 
Fock result (see [Q), which does not correspond to the 
ground states of the systems H + H or H + p. 

Some different interpretation of our result could be ob- 
tained when we consider two separate hydrogen atoms 
with anti-parallel spins of electrons and trace their wave 
functions when shortening the distance between them. 
Let us examine, for example, a dipole moment (z at ) of a 
hydrogen atom, when it nears to another hydrogen atom. 
This dependence is presented in Figure |[ As well as the 
overlap integral this value increases nearly exponential 
with reducing R at large distances, but its behavior near 
the critical point is even more pronounced because its 
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FIG. 2: Overlap integral for the spatial parts of the single- 
electron wave functions of the HF hydrogen molecule as a 
function of the internuclear distance (a.u.). 




FIG. 3: Dipole moment of a single electron wave function 
(a.u.) related to the position of the corresponding nucleus. 
The dependence on the internuclear distance. 



derivative becomes infinite here. For R < R CI an evident 
for delocalized electrons relation z at = R/2 takes place. 

It is natural to ask a question about behavior of 
Hartree-Fock wave functions of other states of the hy- 
drogen molecule at large R. A detailed investigation lies 
outside the scope of this communication and we outline 
here the main characteristics of states, which could be 
considered from some points of view as resembling the 
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ground state configurations of H + H or H~ + p in the 
limit R — > oo. In the notation of a united atom they 
are configurations Is f 2po I, Is | 2p f and 2pp, along 
with the ground state Is 2 considered above. There is 
no quantum phase transitions in both singlet and triplet 
configurations ls2p n . As R — > oo the energy of the single- 
determinant wave function for the Is f 2po I (as well as 
for the 2p\ when ipi(p, z) — —ipi(p, —z) and ipi = 1P2) has 
a non-physical limit coinciding with that for Is 2 under 
condition ipi = ifa. On the other hand, due to exchange 
terms the energy of the triplet ls2po configuration has 
the limit corresponding to the ground state of two sepa- 
rate hydrogen atoms with parallel electron spins. A spin- 
symmetrized two-determinant wave function of the sin- 
glet ls2po tends to that of the system H +p as R — > 00. 
The most complicated is the behavior of the configuration 
2j>q, when we do not impose additional symmetries for it. 
It retains its initial atomic symmetries up to R w 2.97. 
For larger R values the wave functions lose their parity 
with respect to the plane z = 0, but they remain equal 
up to R ss 5.9. At this point the second phase transition, 
leading to ipi ^ -02 takes place. However, this series of 
transformations does not lead to a reasonable physical 
result as R — > 00. 

The effect presented above can be considered in a va- 
riety of aspects, and analogies with several circles of phe- 
nomena could be established. One of these items is a 
problem of localization of electrons in molecules (discus- 
sions and references see in ref. [Q, some more recent ap- 
proaches see e.g. in rcf Jl^]). Our result gives a natural 
boundary, separating two regimes: 1. Dclocalized molec- 
ular orbitals and 2. Electrons localized on separate atoms 
(where, e.g. the value 1 — (ipi\i/>2) could be employed as 
a proper characteristic of the localization) . Thus, we can 
consider our result as a solution of this problem for the 
hydrogen molecule, and a possibility of similar solutions 
for more complicated molecules can be expected. 

Another evident analogy are the Mott transition, 
Wigner localization and related topics like the Hubbard 
model. It is possible to consider the effect presented 



above as an analog of the Wigner localization or Mott 
transition on the level of a single molecule. Following 
steps in this direction have to be, of course, investigations 
of systems with more than two electrons, correlated sys- 
tems and so on. Our calculations are carried out in the 
Hartree-Fock approximation which can be considered as 
obsolete and not precise for the hydrogen molecule. On 
the other hand, precise solutions of Schrodinger equations 
for more or less complex systems are unavailable and the 
methods of the Hartree-Fock level are the best of existing 
tools for them. Our result shows that for a simplest non- 
ionized molecule the symmetry breaking in a solution of 
the HF equations allows obtaining physically correct re- 
sults for arbitrary internuclear distances and could give 
a simple insight into various more complicated problems. 
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